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ABSTRACT 

Closed form, approximate functions for estimating the variances and degrees -of- 
freedom associated with the slow crack growth parameters n, D, B, and A * as measured 
using constant stress rate (“dynamic fatigue”) testing were derived by using propagation 
of errors. Estimates made with the resulting functions and slow crack growth data for a 
sapphire window were compared to the results of Monte Carlo simulations. 

The functions for estimation of the variances of the parameters were derived 
both with and without logarithmic transformation of the initial slow crack growth 
equations. The transformation was performed to make the functions both more linear and 
more normal. 

Comparison of the Monte Carlo results and the closed form expressions derived 
with propagation of errors indicated that linearization is not required for good estimates 
of the variances of parameters n and D by the propagation of errors method. However, 
good estimates variances of the parameters B and A ’ could only be made when the 
starting slow crack growth equation was transformed and the coefficients of variation of 
the input parameters were not too large. This was partially a result of the skewed 
distributions of B and A . Parametric variation of the input parameters was used to 
determine an acceptable range for using closed form approximate equations derived from 
propagation of errors. 
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1. INTRODUCTION 


Slow crack growth (SCG) parameters for glasses and ceramics are determined 
by either strength -based or fracture mechanics based test methods. Strength-based 
methods employ smooth test specimens, such as flexural beams or tensile specimens, and 
estimate SCG material parameters from strengths measured over different time intervals. 
Loading is generally done in a static fashion (i.e., “static fatigue”) or in a continuously 
increasing fashion (i.e., “dynamic fatigue”). The strength-based methods are practical 
because the tests are simple, inexpensive, and usually accomplished quickly. 

Strength-based methods directly sample the preexisting flaw distribution within 
or on the surface of the test specimens. Thus the cracks develop from at least some of the 
same sources that are expected to cause failure in 3 component manufactured in a similar 
fashion from the same material. Only the strength-based approaches have been 
standardized [1-5], and, as a result of the critical nature of flight hardware, data for design 
of such components is generally generated with standardized test methodologies. 

The disadvantage of strength-based methods is that the SCG results are subject 
to the scatter inherent in the strength distribution of the material. Thus the estimation of 
SCG parameters from strength data can result in poor statistical reproducibility, and an 
estimate of the parameter variances is very necessary to the design process. In this 
respect, fracture mechanics based approaches, which usually exhibit lower scatter, in 
combination with strength-based approaches might yield the most confidence. 

Closed form, approximate standard deviation functions for the SCG parameter B 
were previously derived for both the “static” and the “dynamic” loading cases [6-9], and 
the accuracy of the solutions were confirmed to some extent via Monte Carlo simulations 
[8]. Varying degrees of success have been reported in using the approximate solutions 
for the design of flight hardware [10, 11]. 

The variance equations derived previously are very general and were determined 
prior to standardization of the associated test methodologies. Thus the functions are not 
particularly convenient for making rapid estimates from the SCG parameters derived with 
current testing standards [3]. Therefore, convenient closed form expressions were 
derived in terms of the SCG parameters determined with one of the standard methods [3] 
and the resulting functions were compared to Monte Carlo simulations. In addition to 
deriving a variance function for the logarithm of the parameter B as done previously [7], 
functions for the variance of the parameter A* and its degree s-of-ffeedom were derived. 
The functions are needed for input to the computer code FLAGRO [12]. 

2. DERIVATION OF CLOSED FORM FUNCTIONS 

One straightforward method for determining the standard deviation (i.e., the 
square root of the variance) associated with a dependant variable is the generation of 
system moments or the law of propagation of errors (POE). It is based on a Taylor series 
expansion of the dependant variable about the means of the independent variables. A 
specific knowledge of the component distributions is not necessary for application of the 
technique. The expansion is usually truncated at the first term, and thus the functions 
being analyzed need to be relatively linear and the CV's (coefficients of variation; i.e., the 
standard deviation divided by the mean) of the independent variables should not be too 
large. A rule-of-thumb that the CV‘s not exceed ~10% has been recommended [13]. 
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Such a CV is relatively small for strength distributions of many ceramics, and 
corresponds to a 2-parameter Weibull modulus of ~12. For nonlinear functions, the 
results can be improved by transforming the function to a more linear space. This results 
in standards deviation parameters in the transformation space (e.g., standard deviation of 
the logarithm of 5, SD !nS , rather than the standard deviation of B, SD B ). In addition, for 
functions with small first derivatives in the range of interest (e.g. large values of n), 
higher order terms in the expansion cannot readily be ignored [14] and POE may become 
inaccurate. If the input variables are normal, then POE estimates are approximately 
normal [14). Thus, any transformation should simultaneously accomplish two goals: 
more linear functions in specific independent variables, and more normally distributed 
independent and dependant variables. The law of propagation of errors implies that the 
estimated value and the associated standard deviation of a system can be estimated from 
[15]: 
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where E(x l ) is the expected value of variable x t (i.e. the mean), SD Xi is the standard 
deviation of variable x,-, and Covfaj , Xj) is the covariance between x, and xj. Any number 
of higher order terms can be included in estimates of the standard deviations; however, 
estimation of the moments associated with them is cumbersome, and terms greater than 
second order in equation (2) are frequently dropped for convenience. Also, the second 
term in equation (1) is usually dropped. This makes transformation an appealing 
approach for improving the accuracy of the technique. 

2.1 Constant Stress Rate Testing 

For most ceramics and glasses, the slow crack growth rate above the slow' crack 
growth limit is expressed by the following power-law relation: 

V -~Z~ AK" =A* [Elj" (3) 

dt Kic 

where v, a and t are crack velocity, crack size and time, respectively. A and n are the 
material/environment dependent SCG parameters and Kj and K IC are, respectively, the 
Mode I stress intensity factor and the critical stress intensity factor or fracture toughness of 
the material. 

For constant stress rate or “dynamic fatigue” testing, the corresponding fracture or 
fatigue strength. Of, is expressed as a function of stress rate, & , as follows [16]: 
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Gf = [fl(« + i)cr; 2 a\ I/{n+l} (4) 

where O is the expected value of the inert strength and B is a parameter associated with 
A, n, fracture toughness, and the geometry correction factor, Y, for the stress intensity 
factor: 
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The American Society for Testing and Materials (ASTM) approved a full- 
consensus standard [3] for estimating the parameters n and D from equations (4) and (6). 
However, several design codes [12, 17] require not only the parameter n, but either the 
fatigue parameter B or A as defined in equations (5) and (3). Further, as mentioned 
previously, the standard deviations of these parameters are necessary in estimating 
confidence intervals on SCG predictions and it is convenient to calculate the parameters and 
standard deviations directly from the regression statistics of the linearized form of equation 
(4) as defined by the ASTM standard [3]. 

In the ASTM standard, estimates of the SCG parameters are determined from 
the slope, a. and the intercept, ft, of a plot of log ]0 a f versus log L o & by writing 
equation (4) as 

log 10 G f = — tog 10 <7 + log l0 D ( 6 ) 

n + 1 

where 

log lri D = - I -log l0 [B(n + l)ar 3 ]. (7) 

n + 1 

In the linear regression analysis, the dependent variable is logio Of and the independent 
variable is logio & . 

2.2 Parameters Determined Without Transformation 


For constant stress-rate or “dynamic fatigue” testing, application of equations (1) 
and (2) without transformation gives the resulting SCG parameters: 


Estimated n : 

Estimated Standard Deviation of n : 


n = 1 

a 

(8) 

SD„ 
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Estimated D: 
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Estimated Standard Deviation of D: 


SD d » (in 1 0)(SDp )(l0^ ) (11) 
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Estimated B: 


( 12 ) 



where Q = a - ft In 10 + In (S i . 

The statistics for parameter^’ corresponding to equation (3) can be calculated from 



The use of equation (2) for deriving standard deviation functions assumes that a single 
term in the Taylor series expansion is sufficient. In addition, it was assumed that only the 
regression parameters a and /? are correlated, and that the estimated parameters can be 
substituted for the expected values. The “hats” and “bars” used to indicate “estimate" 
and “mean” (e.g., p and CT) have not been included for brevity. 

The derivations of equations (12) and (14) only consider the first term in 
equation (1). If all three terms of (1) are retained, the following functions are obtained: 
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23 Parameters Determined With Transformation 


Equation (1) becomes exact in the case of the dependent variable being a linear 
function of the independent variables. Equations (12) and (14) were transformed by 
equating the logarithms of each side: 



For the purposes of differentiation, the equations were taken as linear in a, p lnKj c , and 
ln<Ji. If all the terms in equation (1) are used to estimate the parameters, then 
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Application of equation (2) to (1 8) and (1 9) gives 
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Tile estimated standard deviations of InA* and InB differ from the estimated standard 
deviations of A* and B nominally by factors of A * and B, respectively. 


3. Monte Carlo Simulation 

The computing power available today makes Monte Carlo simulation an 
appealing technique for parameter estimation. In this technique, values for the 
independent variables are chosen at random from specified distributions and a value of 
the dependant parameter is calculated. Thus, knowledge of, or assumptions about the 
distributions of the independent parameters is required. The technique is repeated a large 
number of times, and the statistics of the resultant data pool are calculated. The 
simulations in this work were repeated 50,000 times in order to insure a reasonable 
sample size at common probability limits. 

3.1 Input Data 

Data for input to the simulation came from fracture toughness, inert strength, 
and constant stress rate SCG testing of sapphire for a window application [18]. Linear 
regression of the data according to equation (6) resulted in the parameters in Table 1. 
The distributions of a, /? and K Jc were assumed to be normal. Although the distribution 
of O' is treated as normal, values less than zero are physically unrealistic as they imply no 
stress corrosion, and the distribution could be truncated for values less than zero. 
Additionally, values of greater than -0.1 imply that the assumption of n > -10 used in 
derivation of the slow crack growth laws [16] is violated. To determine the significance 
of truncation, Monte Carlo runs were performed with and without truncation of the 
distribution of or below 0.01 and below zero. No differences were noted in the estimates 
and standard deviations. Truncation was not performed in subsequent runs. In addition 
to the SCG and fracture toughness data, inert strength data was generated by testing in 
high purity dry nitrogen. The strength data was fit to both 2-parameter and 3-parameter 
Weibull distributions. The resultant parameters are summarized in Table 2. 
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Table 1. Estimated means and standard deviations used for normal distributions. 


Parameter 

Mean 

Standard 

Deviation 

Coefficient of 
Variation, % 

Estimated slope, a 

0.05341 

0.00367 

6.9% 

Estimated intercept, ft 

2.9131 

0.00654 

0.2% 

Fracture toughness, K ja MPaVm 

2.51 

0.12 

4.8% 

Geometry Correction, Y 

1.951 

0 

0 

Inert Strength, <r h MPa 

1108 

93.5 

8.4% 


Table 2. Weibull parameters for inert strength distributions. 


Inert Strength 
Weibull Parameters 

Threshold 
Strength, MPa 
<y,h 

Characteristic 
Strength, MPa 

Weibull 

Modulus 

m 

Two Parameter 

0 

1148.8 

14.6 

Three Parameter 

738.2 

359.0 

3.89 


3.2 Estimation of Covariance 

The dependence of the value of ft upon the value for a was included in the 
Monte Carlo simulations by using a closed-form covariance term. Many existing 
computer codes calculate regression statistics such as the standard error of the parameters 
(i.e. the square root of the variances). Thus the covariance, Cov(a , f$), can easily be 
calculated from 

Cov{a, fi)=-V {ctyog a (24) 


where V(a) is the variance of the regression slope (i.e. the standard error squared), and 
log o- is the mean of the log of the stressing rates applied. Note that the ASTM standard 
[3] uses base 10 logarithms for plotting and regression, and thus equation (24) uses base 
1 0. The covariance between the two can also be found from [7,8]: 
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and a is the standard error of regression, & y is the f h value of stressing rate, and J is the 
number of data points. 


4. RESULTS 

4.1 Results without Transformation 

The distributions for n, D, B, and Al from the Monte Carlo simulation are shown 
in Figures 1-4. The distributions for n and D are relatively symmetric, whereas the 
distributions for B and A* are extremely skewed to the left. 

The estimates of the parameters n, D, B, and A* and their variances from the 
Monte Carlo simulation are compared to the functions that were derived with the law of 
propagation of errors in Tables 3-6. This comparison shows good agreement betw een the 
closed-form estimates and the Monte Carlo results for the parameters n and D, as 
summarized in Table 3. 


Table 3. Comparison of estimates for parameters n and D 


Parameter 

Monte Carlo 

Monte Carlo 
(repeat run) 

POE Estimate 

Wmean 

17.81 

17.81 

17.72 

ft median 

17.76 

17.72 

SD„ 

1.285 

1.306 

1.286 

&mean 

819.2 

818.7 

818.7 

& median 

818.9 

818.5 

SD d 

12.34 

12.32 

12.34 
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n 

Figure 1: Histogram for slow crack growth parameter n from Monte Carlo simulation based on constant stress 
rate testing of sapphire in water. 
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D 

Figure 2: Histogram for slow crack growth parameter D from Monte Carlo simulation based on constant stress 
rate testing of sapphire in water. 
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Figure 4: Histogram for slow crack growth parameter^* from Monte Carlo simulation based on constant stress 


rate testing of sapphire in water. 
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For B and A*, the POE solutions and the Monte Carlo simulation agree only when median 
values are calculated from the Monte Carlo data. This is a result of the distributions of A 
and B being very skewed. The skewness is the result of the highly nonlinear nature of the 
functions and the associated sensitivity to small variations in the material properties, 
especially the inert strength. Note that the median values from the 3-parameter Monte 
Carlo match well with the POE solutions when the mean inert strength is used in the 
closed form solutions, and that the median values of the 2-parameter Monte Carlo match 
well with the POE solutions when the median inert strength is used with the POE 
solutions. This occurs because the 3-parameter Weibull distribution is relatively 
symmetric, whereas the two-parameter Weibull distribution exhibits a large shape 
parameter and a more skewed distribution. In general, the standard deviations from the 
two approaches do not agree. 


Table 4. Comparison of POE and Monte Carlo estimates for parameter B. 


Parameter 

POE 

(mean c-) 1 

POE 

(median oj) 2 

Monte Carlo 
3-parameter 
Weibull 

Monte Carlo 
2 -parameter 
Weibull 

Parameter via eq. (12); SD via eq. (13) without the Cov term: 

fttnean 

251,203 


742,130 

251,011 

ft median 

— 

219,009 

248,387 

204,666 

SD b 

359,656 

312,025 

1,695,835 

20,138,202 

Parameter via eq. (12) plus a Cov term ; SD via eq. (1 3): 

ft mean 

268,650 


807,512 

1,872,491 

ftmedian 


234,651 

264,832 

236,521 

SD b 

373,709 

324,638 

1,838,818 

32,624,604 

Parameters with both second order terns (eq 16): 

ftmean 

513,825 


1,730,407 

5,709,024 

ftmedian 


444,878 

512,423 

443,339 


1 . Mean value of a; = 1 1 08.05 MPa. 

2. Median value of o; = 1 1 17.8 MPa 
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Table 5. Comparison POE and Monte Carlo results for parameter A *. 


Parameter 

POE 
(mean o;) 1 

POE 

(median o; ) 2 

Monte Carlo 
3-parameter 
Weibull 

Monte Carlo 
2-parameter 
Weibull 

Parameter via eq. (14); SD via eg. (15) without the Cov term: 

a' 

si mean 

8.381E-07 

— 

2.175E-06 

1.863E-06 

-A median 

— 

9.613E-07 

8.494E-07 

9.618E-07 

SDA* 

1.1836E-06 

1.349E-06 

4.721E-06 

3.1 12E-06 

Parameter via eq. (14) plus a Cov term; SD via eq. (15): 


9.047E-07 


2.782E-06 

2.071E-06 



1.040E-06 

9.223E-07 

1.042E-06 

SD a * 

1.222E-06 

1.395E-06 

8.180E-06 

3.644E-06 

Parameters with both second order terms (eq. 17): 

A mean 

1.712E-06 


4.662E-06 

4.003E-06 

A median 


1 .957E-06 

1.707E-06 

1.947E-06 


1. Mean value of <7j - 1 108.05 MPa. 2. Median value of tJ, = 1 1 17.8 MPa 


4.2 Results with Transformation 

The results in Tables 4 and 5 indicate that POE estimates for B and A* are not 
acceptable without transformation. The resultant distributions for InB and InA* are shown 
in Figures 5 and 6, and the results from POE and Monte Carlo simulations are compared 
in Tables 6 and 7. Logarithmic transformation of the SCG equation brings the two 
methods into good agreement. As seen earlier, the Monte Carlo simulation using a three- 
parameter Weibull distribution compares well with POE results that use the mean value 
of the inert strength. And the Monte Carlo simulation using a two-parameter Weibull 
distribution matches the closed form results calculated with the median inert strength. 
The addition of the second order terms in equation (1) lowered the estimated B and A by 
~4%. Generally the standard deviations from POE are ~5% smaller than those from 
Monte Carlo. 
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Figure 5: Histogram for slow crack growth parameter tn(B) from Monte Carlo simulation based on constant 
stress rate testing of sapphire in water. 
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Figure 6: Histogram for slow crack growth parameter In(A’) from Monte Carlo simulation based on constant 
stress rate testing of sapphire in water. 
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Table 6. POE and Monte Carlo results for parameter lrt(B). 


Parameter 

POE 

(mean o}) 1 

POE 

(median oj) 2 

Monte Carlo 
3 -parameter 
Weibull 

Monte Carlo 
2-parameter 
Weibull 

Parameter via eg. (18); SD via eq. (22) without the Cov term: 

Ifl (13) mean 

12.43 

— 

12.46 

12.45 

If! (B) median 

— 

12.30 

12.44 

12.32 

SDln(B) 

1.452 

1.456 

1.518 

1.553 

Parameter via eg 

. (18) plus a Cov term; SD via eq. (22): 

■ "_'2 

12.42 


12.43 

12.43 

I” (B) median 


12.28 

12.41 

12.30 

msmm 

1.507 

1.512 

1.511 

1.563 

Parameters with both second order terms (eq. 20): 

ln(B) mecm 

12.39 

— 

12.43 

12.43 

In (B) median 


12.25 

12.41 

12.30 


1 . Mean value of <7 t = 1 1 08.05 MPa 2. Median value of a t = 1 1 17.8 MPa 


Table 7. POE and Monte Carlo results for parameter In (A*) 


Parameter 

POE 

(mean cj; ) 1 

POE 

(median <7i ) 2 

Monte Carlo 
3 -parameter 
Weibull 

Monte Carlo 
2-parameter 
Weibull 

Parameter via eq. (19); SD via eq. (23) without the Cov term: 

ln(A )mean 

-13.99 


-14.02 

-14.01 



-13.85 

-14.00 

-13.87 

SDi„(a *) 

1.449 

1.451 

1.488 

1.524 

Parameter via eq. 

f 19) plus a Cov term; SD via eq. (23): 

In (A )mean 

-13.98 


-13.99 

-13.99 



-13.84 

-13.97 

-13.85 


1.478 

1.482 

1.480 

1.534 

Parameters with both second order terms (eq. 21): 


-13.95 

— 

-13.99 

-13.99 

In (A ) median 


-13.81 

-13.97 

-13.85 


1 . Mean value of o; = 1 1 08.05 MPa 2. Median value of c; - 11 1 7.8 MPa 
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It is noteworthy that the value of B estimated in log space (i.e. for lnB= 12.43, B = 
250,196) either with or without the second order terms of equation (20) is very similar to 
the direct estimate of B (251,203) only when the second order terms are not included in 
the direct estimate. If the second order terms are included, the direct estimate of B 
becomes twice as large (251,203 vs. 513,825). Evidently, the best estimates of B and A' 
without linearization are generated by using the initial algebraic formula and Cov. 

4.3 Generation of Confidence Bands 

For a small data set (< 40), probability limits are placed on parameters by using 
the Student’s t distribution: 

B Upper = EXP[lnB ± t{SD , na )] . (27) 

Lower 


Calculation of t requires the DOF (degrees of freedom) <p, which is not readily estimated 
for a complex system. Welch [19] gave a formula for estimating the DOF when multiple 
components of variation are involved in comparison of two values. Jacobs and Ritter [8] 
generalized the Welch formula for estimation of the DOF of populations with different 
variances: 



(28) 


where tps represents the DOF for the system variable of interest, SD, represents the 
component of variation contributed to the system by variable i, and <p, represents the DOF 
of the components adding variation to the system. For SCG parameters, the two 
components adding variation are the inert strength and the regression parameters a and (3. 
Application of (28) to (22) leads to 
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and results in 47 DOF in SD lnB base on 14 data points in inert strength (0* = 13) and 41 
data points in the regression (0^ = 39). Application of equation (28) to equation (23) 
leads to 
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If, however, sufficiently large data sets are used, then the normal distribution can 
reasonably be assumed, and two standard deviations correspond to 95% confidence, etc: 

^Upper95% — EXP[ln B ± 2SD, rI1 ] (31) 

Lawer9S% 


Probability limits calculated on the parameters with equation (29) and from the Monte 
Carlo data are summarized in Tables 8 and 9. The POE and Monte Carlo 3 -parameter 
results agree reasonably. However, the 2-parameter Monte Carlo results are less 
conservative (failure time is proportional to B) and result in larger values for B and 
smaller values for A*. This appears to result from a more skewed distribution of InB that 
results in larger 95% values. Use of the median strength in the POE gives the most 
conservative result (i.e. the smallest B and largest A' values). It is noteworthy that the 
standards for determining the inert strength and the SCG parameters recommend using 30 
and 40 test specimens, respectively. Thus equation (31) is likely sufficient for most 
purposes. As the least value of B will produce the shortest life, a one-sided limit may be 
more appropriate in design than the two-sided limits presented in Tables 8 and 9. 


Table 8. POE and Monte Carlo 95% probability limits for parameter B based on 
equations (1 8), (20), (22) and (27). 



POE 

(mean a) 1 

POE 

(median cJj ) 2 

Monte Carlo 
3 -parameter 
Weibull 

Monte Carlo 
2-parameter 
Weibull 

Parameters via eq. (18); SD via eg. (22) without the Cov term: 

Lower 

Bound 

12,258 

10,601 

13,202 

16,448 

Upper 

Bound 

5,147,88 5 

4,519,110 

5,003,181 

8,361,753 

Parameters with second order terms (eq. 20); SD via eq. (22): 

Lower 

Bound 

10,485 

9,038 

14,885 

18,528 

Upper 

Bound 

5,540,317 

4,871,973 

5,1278,769 

8,509,182 


1 . Mean value of a t = 1 108.05 MPa. 2. Median value of 0/ = 11 17.8 MPa 
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Table 9. POE and Monte Carlo 95% probability limits for parameter A * based on 
equations (19), (21), (23) and (27). 



POE 

(mean o; )’ 

POE 

(median 0 ; ) 2 

Monte Carlo 
3 -parameter 
Weibull 

Monte Carlo 
2-parameter 
Weibull 

Parameters via eq. (19); SD via eq. (23) without the Cov term: 

Lower 

Bound 

3.750 x 10~ 8 

4.276 x 1 0' 8 

4.286 x 10’ 8 

2.567 x 10‘ 8 

Upper 

Bound 

1.873 x 10’ 5 

2.160 x 10‘ 5 

1.459 x 10~ 5 

1.163 x 10’ 5 

Parameters with second order terms (eq. 21); SD via eq. (23): 

Lower 

Bound 

3.661 x 10' 8 

4.166 xlO' 8 

4.438 x 10‘ 8 

2.637 x 10' 8 

Upper 

Bound 

2.075 xlO 5 

2.405 x lO’ 5 

1.665 x 10’ 5 

1.328 x 10 5 


1 . Mean value of o; = 1 1 08.05 MPa. 2. Median value of Oi = 1 1 1 7.8 MPa 


4.4 Range of Applicability for POE 

To determine the effect of input variance on the accuracy of the equations 
derived from POE, the CV’s of a, P. and <7, were varied above that observed, and the 
standard deviation of InB was calculated. In order to calculate the SDj nB , the standard 
deviation of InCTi is needed. This was calculated from a Monte Carlo simulation that used 
the characteristic strength and a 2-parameter Weibull modulus corresponding to the 
desired CP as calculated with the approximation m ~ 1.2/CP [7, 20]. The strength data 
resulting from the Monte Carlo simulation was converted to InOi and the corresponding 
SDi„ a was calculated. 

The ratio of the standard deviations from the POE and the Monte Carlo 
approaches are compared in Table 10. Generally, the POE method produced smaller 
values that were within -1 5% of those from Monte Carlo. As might be expected, better 
comparisons occur for small CV’s of a and p. Surprisingly, the best results occur for 
large and intermediate standard deviations of the inert strength. This may be the result of 
the logarithmic transform being excessive. 

To investigate the applicability of the equations for larger n values, the 
sensitivity for an alumina with n - 47 [21] was determined. The distribution or of was 
truncated for values of a < 0.008. Again, relatively good comparisons were obtained. 
However, it is noteworthy that the material had a large Weibull modulus (-“25) and low 
CV a and CVp. Despite the small initial CP’s, the size of the estimated SD !nB relative to 
InB tends to become quite large as the input CP for the inert strength is increased to 
values common for ceramics, and the mean and median values of InB as calculated via 
Monte Carlo tend to diverge. For example, with CV a = 15%, CV p - 0.5%, and CV a ~ 
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30%, equations (18) and (20) give InB - -5.67 and InB = -6.14 respectively with the mean 
strength, and InB - -5.93 and InB = -6.41 respectively with the median strength. Whereas 
Monte Carlo simulation is much less conservative, giving a mean value of InB - -0.46 
and a median InB = -2.21. In this case, despite the transformation to logarithmic space, 
the distributions of InB and InA are skewed as compared to those for n — 18. More 
investigation with large values of n is needed. 


Table 10. Ratio of 5Z>/^ as calculated from POE via equation (22) and Monte Carlo 
simulation with m ~ 1.2/CF. 


SDg(poE)! SD b<mc) 

cv a = 

10% 

10% 

20% 

20% 

CV <T t -* ~ m 

CVp= \ 

0.5% 

1.0% 

0.5% 

1.0% 

10%-* m « 12 

0.94 

0.92 

0.82 

0.77 

15% 8 

0.96 

0.94 

0.87 

0.82 

20% ^ m ~ 6 

1.00 

0.98 

0.91 

0.87 

30% “*• 4 

0.98 

0.97 

0.91 

0.89 


Table 11. Ratio of SD /nB as calculated from POE via equation (22) and Monte Carlo 
simulation with m ~ 1.2/CF. Calculations are for an alumina [21] with n - 47.3, a- 
0.0207 ± 0.0020, 0- 2.2747 ± 0.0048, a median inert strength of 278.6 MPa, and a mean 
inert strength of <% = 277.0 MPa ± 12. 


SD B( poe)I SD B (mq 

cv a = 

10% 

10% 

15% 

15% 

CV a . —*■ ~m 

CVp = 

0.25% 

0.5% 

0.25% 

0.5% 

\0%~* m - 12 

0.99 

0.98 i 

0.95 

0.92 

15% “*** m - 8 

1.00 

0.99 

0.96 

0.95 

20% m ~ 6 

1.00 

0.98 

0.97 

0.95 

30% — ► m~4 

0.98 

0.98 

0.96 

0.96 
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5. CONCLUSIONS 


Two methods were used to approximate variances of SCG data: Propagation of 
errors and Monte Carlo simulation. Transformation of the functions describing the 
parameters is necessary to reasonably predict the standard deviations of parameters InB 
and In A*. The addition of a covariance term slightly improved estimations of the 
standard deviations of ln(B) and In (A*) for two cases considered (n - 18 and n - 47). 

The addition of second order terms made little difference for POE estimation of 
the parameters ln(B) and ln(A*) for sapphire with n = 18. For an alumina with n = 47, the 
differences were significant, and the POE estimates of InB and InA* tended to be very 
conservative relative to Monte Carlo simulations. 

Propagation of errors is a useful approximation for estimating variances of slow 
crack growth parameters, however, it is an approximation that results in estimates within 
-15% of Monte Carlos simulations for the range of variances commonly encountered in 
SCG testing of ceramics. Further study with large n values is needed. Equations (8), 
(11), (22), and (23) provide reasonably accurate functions for determining standard 
deviations of slow crack growth parameters determined from constants stress rate testing. 
For very critical applications, it may be more appropriate to use Monte Carlo or 
Bootstrap simulations, especially considering the availability of computing power. 
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